Effect of the magnetic resonance on the electronic spectra of high Tc superconductors 
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We explain recent experimental results on the superconducting state spectral function as obtained 
by angle resolved photoemission, as well as by tunneling, in high Tc cuprates. In our model, 
electrons are coupled to the resonant spin fluctuation mode observed in inelastic neutron scattering 
experiments, as well as to a gapped continuum. We show that, although the weight of the resonance 
is small, its effect on the electron self energy is large, and can explain various dispersion anomalies 
seen in the data. In agreement with experiment, we find that these effects are a strong function of 
doping. We contrast our results to those expected for electrons coupled to phonons. 
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I. INTRODUCTION 

Understanding superconductivity in the cuprates is 
one of the great challenges of physics. Determining 
the nature of single particle excitations is of fundamen- 
tal importance for achieving this goal. Two types of 
experiments have been extensively used to study such 
excitations: angle resolved photoemission spectroscopy 
(ARPES) and tunnehng. 

In this paper, which deals with the superconducting 
state only, we address the questions, what the spectral 
properties of fermionic excitations are, and how their low- 
energy dispersion is renormalized. We do not directly ad- 
dress the question of the origin of superconductivity in 
the cuprates. Rather, we assume that an effective pair- 
ing interaction exists, and study the additional effects 
which coupling to certain collective excitations present 
in cuprates have in renormalizing single particle proper- 
ties. The corresponding collective excitations responsible 
for such renormalizations are most directly seen in other 
types of experiments. One of them, inelastic neutron 
scattering, gives the most useful information about both 
phonons and magnetic excitations in the energy range of 
interest (< lOOmeV). 

Motivated by earlier workipiSiiii^iSJ we have presented 
in Ref.^a model which describes the ARPES and tunnel- 
ing spectra. Here, we describe details of our calculations, 
and extend them by including the effect of the spin fluc- 
tuation continuum. In addition, we address the issue of 
the doping dependence of the ARPES spectra. Finally, 
for comparison, we discuss the effect on the electrons of 
coupling to a particular phonon which was recently sug- 
gested to account for the renormalization of the ARPES 
dispersion in the nodal regions of the zone. 

Our outline is the following: starting in Section II from 
the information which experiments give about single par- 
ticle properties of low lying excitations in cuprates, we 
look for a suitable collective excitation which best fits 
the data. Then, we develop in Section III a model in 
which the collective mode is identified as the magnetic 
resonance observed in inelastic neutron scattering exper- 
iments. The results of calculations using this model are 
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FIG. 1: Notation used for special points in the Brillouin zone. 
The Fermi surface is shown as a black curve. Equal energy 
contours are shown in gray for energies between ±50meV. 
The dispersion used here was obtained by a 6-parameter tight 
binding fit to angle resolved photoemission dispersions in opti- 
mally doped Bi2Sr2CaCu208_5.* The dispersion has a saddle 
point at the M point. The A*" point corresponds to the node 
of the d-wave order parameter in the superconducting state. 



presented in great detail. Finally, in Section IV, we ad- 
dress the question what electron-phonon coupling con- 
tributes to renormalization effects on the dispersion. Sec- 
tion V offers a brief summary. 



II. EXPERIMENTAL EVIDENCE 
A. Angle resolved photoemission 

It has been known for some time that near the 
(tt, 0) (M) point of the zone, the spectral function in 
the superconducting state of Bi2Sr2CaCu208+5 shows 
an anomalous lineshape, the so called 'peak-dip-hump' 
structure^iSiifliii This structure was also found recently 
in YBa2Cu307-5r^ and in Bi2Sr2Ca2Cu30io-H5iii^ For 
the notation of special points in the Brillouin zone which 
we use throughout this paper, see Fig. ^ 

Extensive studies on Bi2Sr2CaCu208+5 as a function 
of temperature revealed that this characteristic shape of 
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the spectral function is closely related to the supercon- 
ducting state. In the normal state, the ARPES spec- 
tral function is broadened strongly in energy, the broad- 
ening increasing with underdoping.™ The width of the 
spectral peak quickly decreases with decreasing temper- 
ature below and sharp quasiparticle peaks were 
identified well below Tc along the entire Fermi surface.^® 
When lowering the temperature below Tc, the coher- 
ent quasiparticle peak grows at the position of the lead- 
ing edge gap, and the incoherent spectral weight is re- 
distributed to higher energy, giving rise to a dip and 
hump structure4*2iifl This peak-dip-hump structure is 
most strongly developed near the M-point of the Bril- 
louin zone. The well defined quasiparticle peaks at low 
energies contrasts to the high energy spectra, which show 
a broad linewidth which grows linearly in energy. 
This implies that a scattering channel present in the 
normal state becomes gapped in the superconducting 
stateJ^ The high energy excitations then stay broadened, 
since they involve scattering events above the thresh- 
old energy. While this explains the existence of sharp 
quasiparticle peaks, a gap in the bosonic spectrum which 
mediates electron interactions leads only to a weak dip- 
like featured This suggests that the dip feature is in- 
stead due to the interaction of electrons with a sharp (in 
energy) bosonic mode. The sharpness implies a strong 
self energy effect at an energy equal to the mode energy 
plus the quasiparticle peak energy, giving rise to a spec- 
tral dip.^ The fact that the effects are strongest at the 
M points implies a mode momentum close to the (7r,7r) 
wavevector^ 

More clues are obtained by studying the dispersion 
of the related self energy effects. Recent advances in 
the momentum resolution of ARPES have led to a de- 
tailed mapping of the spectral lineshape in the high Tc 
superconductor Bi2Sr2CaCu208+i throughout the Bril- 
louin zoneiSii2S The data indicated a seemingly unrelated 
effect near the d-wave node of the superconducting gap, 
where the dispersion shows a characteristic 'kink' fea- 
ture: for binding energies less than the kink energy, the 
spectra exhibit sharp peaks with a weaker dispersion; be- 
yond this, broad peaks with a stronger dispersion iiSi2ii2S 
This kink is present at a particular energy all around the 
Fermi surface?^ and away from the node, the dispersion 
as determined from constant energy spectra (momentum 
distribution curves, MDCs) shows an S-like shape in the 
vicinity of the kink^ The similarity between the exci- 
tation energy where the kink is observed and the dip 
energy at M, however, suggests that these effects are 
related.* Additionally, the observation that the spectral 
width for binding energies greater than the kink energy 
is much broader than that for smaller energiesiSiSiiS^ is 
very similar to the difference in the linewidth between 
the peak and the hump at the M points. Further ex- 
perimental studies supported the idea of a unique energy 
scale involvedfS They found that away from the node, 
the kink in the dispersion as determined from constant 
momentum spectra (energy distribution curves, EDCs) 



develops into a 'break'; the two resulting branches are 
separated by an energy gap, and overlap in momentum 
space. Towards M, the break evolves into a pronounced 
spectral 'dip' separating the almost dispersionless quasi- 
particle branch from the weakly dispersing high energy 
branch (the 'hump'). The kink, break, and dip features 
all occur at roughly the same energy, independent of po- 
sition in the zone,^^ the kink being at a slightly smaller 
energy than the break featurei^ 

The high energy dispersion is renormalized up to at 
least 200 meV and does not extrapolate to the Fermi 
surface crossing i2i42^ This lets us conclude that the con- 
tinuum part of the bosonic spectrum coupling to the 
fcrmionic excitations extends to high energies. 

Finally, there is important information contained in 
the doping dependence of the self energy effects. In un- 
derdoped compounds, there is a pseudogap between Tc 
and the pseudogap is maximal near the M-point 

of the Brillouin zone and is zero at arcs centered at the 
A^-points which increase with temperature^ In the pseu- 
dogap state above Tc, there are low energy renormaliza- 
tions in the dispersion, and some trace of the kink feature 
persists. But in the recent work by Johnson et al.'^^, it 
was clearly shown that an additional renormalization of 
the dispersion sets in just at Tc. This indicates that the 
bosonic spectrum redistributes its spectral weight when 
entering the superconducting state. The additional low 
energy renormalization of the dispersion below the kink 
energy follows an order parameter like behavior as a func- 
tion of temperaturci^ Arguing that the renormalization 
near the nodal regions is influenced by the coupling to the 
same bosonic mode which causes the strong self energy 
effects at the M point of the Brillouin zone, the above 
implies that some mode intensity may be present in the 
pseudogap state already, but there is an abrupt increase 
in the mode intensity when going from the pseudogap 
state into the superconducting state, and this increase 
shows an order parameter like behavior as a function of 
temperature below Tc. 

The energy of the mode, as inferred from the en- 
ergy separation Qq between the peak and the dip, was 
shown to decrease with underdoping.^* Similarly, the 
kink energy is maximal at optimal doping and de- 
creases both with underdoping and overdopingj24 indi- 
cating some relationship between the kink at the nodal 
N point and the peak-dip-hump structure at the M 
point. With underdoping, the sharp quasiparticle peak 
moves to higher binding energy, indicating that the gap 
increases.^* At the same time, the spectral weight z of the 
peak drops^*'^^ leaving the quantity zIS.m /ksTc roughly 
constant. "^"^ Also, the hump moves to higher binding en- 
ergy and loses weight with underdoping.^'* This doping 
evolution of the quasiparticle peak points to an increasing 
mode intensity at the (tt, tt) wavevector with underdop- 
ing. Again, there is a similarity to the nodal direction: 
the low energy renormalization of the dispersion below 
the kink energy increases with underdoping,^^ consistent 
with a common origin of both effects. 
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B. C-axis tunneling spectroscopy 

Unusual spectral dip features in tunneling data of 
Bi2Sr2CaCu208-5 are found in point contact junctionS)2i 
in scanning tunneling spectroscopy (STM) )?^i?? in break 
junctionSj-^''^^ and in intrinsic c-axis tunnel junctions^ 
Consistently, these data show a peak feature, usually as- 
signed to the maximal c?-wave superconducting gap, and 
a hump feature at higher bias, separated from the peak 
by a pronounced dip feature. A characteristic of this dip 
feature in SIN junctions is that it occurs asymmetrically 
around the chemical potential, usually stronger on the 
occupied side of the spectrum i^iiSSiS This asymmetry 
was succesfully explained within the theoretical model 
presented below i& The dip feature has been observed in 
tunneling spectra of the single CU-O2 layer compound 
TlzBaaCuOe as well.^ 

In order to extract information about the bosonic mode 
which would produce a dip feature in the tunneling con- 
ductance, a systematic study as a function of doping was 
performed in break junction tunneling spectroscopy by 
Zasadzinski et al.^"^ There, the doping dependence in 
Bi2Sr2CaCu208_5 of the peak-dip-hump structure was 
determined over a wide range of doping. It was found 
that the dip-peak energy separation, follows Tc as 
Q, = 4.9fcsTc. As expected for an excitonic mode, H. ap- 
proaches but never exceeds 2A in the overdoped region, 
and il/A monotonically decreases as doping decreases 
and the superconducting gap increases. The dip feature 
is found to be strongest near optimal doping. Similar 
shifts of the dip position with overdoping were reported 
prevously by STM.^^ Together with the ARPES results, 
these studies give a detailed picture about the doping 
dependence of the mode energy involved in electron in- 
teractions in the superconducting state. 



III. COUPLING TO THE MAGNETIC 
RESONANCE MODE 

There have been several theoretical treatments which 
assigned the anomalous ARPES lineshape near the M 
point of the zone to the coupling between spin fluc- 
tuations and electrons i2424i^iSi& The 'collective mode 
model' proposed by Norman and Ding^ was suggested 
to account for the unusual APRES lineshapes by cou- 
pling electrons to a dispersionless collective mode. The 
main motivation for a more detailed study of this model 
m Ref. ^ was to additionally account for the dispersion 
anomalies ('kink'), and the isotropy and robustness of 
this characteristic energy scaleiS^ 

The minimal set of characteristics for the collective 
mode we are interested in follows from the experimen- 
tal results from ARPES and tunneling. The mode 
is characterized by its energy and its intensity at the 
(tt, it) wavevector (the wavevector being suggested by the 
momentum dependence of the strength of the ARPES 
anomalies). Its properties from ARPES and SIS tunnel- 



ing are as follows. The energy should be weakly depen- 
dent on momentum, roughly 40 meV in optimally doped 
cuprates, follow Tc with doping, and be constant with 
increasing temperature up to Tc- The intensity should 
be maximal at the (tt, tt) wavevector, where it should in- 
crease with underdoping and follow an order parameter 
like behavior as a function of temperature below Tc- The 
mode should be absent in the normal state; a remnant 
can be present in the pseudogap state, but an abrupt 
increase in intensity should occur at Tc with lowering 
temperature. 

A sharp resonance with characteristics fitting the 
ones extracted from ARPES and tunneling measure- 
ments was observed in inelastic neutron scattering ex- 
periments on bilayer cuprates in the superconduct- 
ing state, with an energy near 40 meV in optimally 
doped compoundsi^2iiMiii2ii^ A similar resonance fea- 
ture at the (tt, tt) wavevector is observed in underdoped 
YBa2Cu306+a:, but at a reduced energy4i^Si^ The 
resonance was also found in Bi2Sr2CaCu208+5, both 
in the optimally doped^MI and overdoped^® regime. 
Recently, the resonance was discovered in single layer 
Tl2Ba2Cu06 compound as well4^ 

To show how well the above criteria fit, we summa- 
rize its characteristics: the resonance is narrow in en- 
ergy and magnetic in origin.^" Its energy width is smaller 
than the instrumental resolution (typically less than 10 
meV) for optimally and moderately underdoped materi- 
als. Strongly underdoped materials show a small broad- 
ening of the order of 10 meV^^^ The resonance lies 
below a gapped continuum, the latter having a signal 
typically a factor of 30 less than the maximum at Q at 
the mode energy.^*^ The mode energy decreases with un- 
derdoping, and has its maximal value of about 40 meV at 
optimal doping4^i^^i^S4i In both underdoped and over- 
doped regimes, the resonance energy, firesi is propor- 
tional to Tc, with Qre. « 5 - S.STc^MLiSi^aia 

An additional aspect, specific to bilayer materials, is 
that it only occurs in the 'odd' channel, which connects 
the bonding combination of the bilayer bands to the an- 
tibonding onei^ The continuum is gapped in both the 
even and odd scattering channels (the even channel is 
gapped by ~60 meV, even in the normal state) iS We 
will address this issue further below. 

The resonance is strongly peaked at the (tt, tt) wavevec- 
tor. The momentum width of the spin fluctuation spec- 
trum is minimal at the resonance energy where it is 
(in contrast to the off-resonant momentum width) only 
weakly doping dependent, with a full width of about 0.22 
j^-ij54j55j56 fpj^jg corresponds to a correlation length $^sfi 
of about two lattice spacings. 

A sharp resonance is not observed above Tc i^i*^ On 
approaching Tc from below, the resonance energy does 
not shift towards lower energy^^i^Siii but its inten- 
sity decreases towards Tc, following an order parame- 
ter like behavior i^242i^24^iS With underdoping, the in- 
tensity at Q = (tt, tt) increases from about 1.6 fi^ for 
YBa2Cu307 to about 2.6 fM^ per unit cell volume for 
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YBa2Cu306.5'^^^ There is clearly an abrupt change 
in resonance intensity at T^, even in underdoped com- 
pounds. 

Note that in underdoped materials, an incommensu- 
rate response develops below the resonance energy^S 
which however never extends to zero energy, but instead 
the spectrum is limited at low energies by the so called 
spin gap Egg.^^ This part of the spectrum behaves differ- 
ently from the resonance part as a function of doping^ 
We will neglect this (weaker) incommensurate part of the 
spectrum in this paper. 

The total spectral weight of the resonance is small 
and amounts to about per formula unit at low 

temperatures4£^ We will show below that the small- 
ness of the weight of the resonance is not an obstacle to 
achieving large self energy effects. 



A. Theoretical model 

We are interested in the renormalizations of the 
fermionic dispersion due to coupling of electrons to a 
sharp spin fluctuation mode at low energies, equal to 
about 40 meV or less. We will assume that supercon- 
ducting order is already established without coupling to 
this resonant feature in the spin flucutation spectrum, 
and thus describe the superconducting state by an in- 
dependent order parameter Afe. This order parameter 
will be chosen to have d-wave symmetry (here and in the 
following the unit of length is the lattice constant a), 



TABLE I: Parameters for the effective dispersion ^fc. 



Afe = AM{coskx - cosky)/2 



(1) 



which takes its maximal value A m at the M point in the 
Brillouin zone. The magnitude Am will be chosen so that 
the calculated peak in the ARPES spectrum at the M 
point, after including self energy effects due to coupling 
to the spin fluctuations, fits the position of the spectral 
peak in experimental ARPES spectra. We stress that 
we do not specify the origin of the pairing interaction re- 
sponsible for the order parameter A^, but the continuum 
part of the spin fluctuations is one of the candidates. We 
also underline that, as our results will show, the spin fluc- 
tuation resonance supports pairing, but does not cause 
superconductivity in and of itself. 

In the model we employ, the retarded Green's func- 
tions, Gf-f., for fermionic excitations in the superconduct- 
ing state is a functional of the normal state electronic dis- 
persion ^fe, the order parameter A^, and the self energies 
due to couphng to spin fluctuations, Sf ^, $f j,. The term 
'normal state' here refers to the state at the same tem- 
perature, but with zero order parameter. We employ a 
six parameter tight binding fit for this dispersion, having 
the form 

cos kx + cos ky 
S.k =tQ+ h h t2 cos k^ COS ky 



+t3 



COS 2kx + COS 2ky 



+ t. 



cos 2kx cos 2ky 



cos 2kx cos ky + cos kx cos 2ky 



(2) 
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-34 meV 


0.8 eV 


0.6 eV 


-0.2376 eV 



Any set of six independent parameters for the dispersion 
determines the parameters to — t^. The six parameters we 
use are the positions of the N (node) and A (antinode) 
points in Fig.Q parameterized by k^N = \kN — ^r| and 
kMA = \kA — ^m|, the band energies at the M and Y 
points, S,M and S,y, the Fermi velocity at the N point, 
vn = I Wat I, and the inverse effective mass along direction 
A/ — F at the AI point, m^/. Table summarizes our 
choices. For reference, the corresponding ti are (eV): 
to = 0.0989, h = -0.5908, = 0.0962, = -0.1306, 
t4 = -0.0507, and = 0.0939. 

The parameter is not known from experiment. We 
set it to a reasonable value to preserve a dispersion shape 
similar to that obtained from band theory. The inverse 
mass at the M point is known to be negative and small in 
the M — T direction, and it was suggested that it could be 
zero, giving rise to an extended van Hove singularity.— 
Here we chose a finite, moderately small value. As we 
show, the inverse effective mass will decrease when cou- 
pling to the spin fluctuation mode is taken into account, 
and it is this renormalized inverse mass which is experi- 
mentally observed. Similarly, the value of the Fermi ve- 
locity at the node is chosen somewhat larger than the 
experimental value, since again, one observes the fully 
renormalized velocity; in our calculation, self energy ef- 
fects renormalize this value to the moderately smaller 
value observed in experiment. When the doping level is 
varied, the band filling varies (to changes), so that the 
van Hove singularity at the M point, ^m, will move rel- 
ative to the chemical potential. Also, the Fermi crossing 
kA moves along the M — Y line. All other band structure 
parameters are expected to be rather insensitive to the 
doping level. 

The 'normal state' dispersion ^fe, and the order pa- 
rameter Afc, are phenomenological quantities, which are 
already renormalized by other effects which we do not 
need to specify, but which are assumed to influence the 
physics only on an energy scale large compared to the 
scale of interest in this paper (50-100 meV). The self en- 
ergies due to spin fluctuations will have a part due to the 
particle-hole continuum, and another part due to the res- 
onance. We will consider two models, a simple form and 
an extended form. In the simple form, we include the 
effect of the continuum part of the spin fluctuation spec- 
trum by a constant renormalization of the normal state 
dispersion and the order parameter. This model will cap- 
ture the main physics for energies below 100 meV, which 
is dominated by the coupling of the electrons to the res- 
onant spin fluctuations. The reason is the following: as 
we will show below, in this energy range the imaginary 
part of the self energies due to the continuum part of the 
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FIG. 2: Left: Self energy for electrons (full lines). The wavy 
line denotes a spin fluctuation. Right: the model spin fluc- 
tuation spectrum we used for the wavy line in the Feynman 
diagram. The mode affects the low energy fermionic proper- 
ties. The continuum part only couples to electrons with higher 
energies, and is neglected in the simple form of the modelf^ 
Damping of electrons at energies above 100 meV is caused by 
the continuum part, and is included in the extended model, 
which we also discuss in this paper. 



spin fluctuations is zero, and the real part (divided by e) 
only varies weakly both in energy and momentum. This 
allows to approximate it by a real constant in that en- 
ergy range, and thus include it into the renormalization 
of ^fc and Afc. For this case, the 'normal state' reference 
is defined as the state with zero order parameter, inter- 
acting with a spin flucutation spectrum having no reso- 
nance part and a continuum part identical to that in the 
superconducting state. The real, physical normal state 
will be different because the spin fluctuation continuum 
changes when going from the normal to the supercon- 
ducting state, leading to an additional renormalization 
of the dispersion. Thus, in the simple form of the model, 
the low energy dispersion which enters the calculations 
will be approximately proportional to the true normal 
state dispersion, but the proportionality factor will not 
be unity. 

At higher energies, the spin fluctuation continuum 
can be excited, and this leads to an additional strong 
fermionic damping. We will study this effect in an ex- 
tended model which explicitly includes the gapped spin 
fluctuation continuum. For this extended model, the 
'normal state' dispersion will have a different renormal- 
ization factor as compared to the simple model above. 
Speciflcally, we use for the extended model the above 
dispersion scaled by the factor 1.5 and shifted back in 
energy, so that stays at its original value of — 34meV. 

We find that all essential features of the self energy 
effects in the superconducting state are obtained using a 
minimal model with a spin fluctuation spectrum shown 
in Fig. 12 

The continuum formally has to be cut-off at high en- 
ergies. This cut-off only affects the real part of the 
self energy, and variation of the cut-off leads to only a 
weakly energy dependent contribution to the renormal- 
ization factor which can be absorbed in the dispersion 
as described above. We discuss the choice of this cut-off 
later. 



The retarded Green's function in spectral representa- 
tion is given as a function of the self energies as. 
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e,k 
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with excitation energies and coherence factors 



e,k 



^e,k 
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2je,.k + |A. 



(5) 
(6) 

(7) 



The renormalized dispersion and gap function are given 
in terms of the diagonal (Sf^) and off-diagonal (<&ffe) in 
particle-hole space self energies, as 
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(8) 
(9) 

= „ "^-"^ (10) 



We will couple electrons to the spin ffuctuation spectrum 
with a coupling constant 5, which we assume to be inde- 
pendent of energy and momentum. The self energies for 
our model are then given in terms of the spectral function 
of the spin fluctuations with energy w and momentum (f, 
Bu>,qi by the expressions (we chose a representation espe- 
cially well suited for numerical studies, see App. 

— \^ ^2 n ^R 

^TJ2GtL,i^e,.)9'D'/ie-^e,,) (11) 



-TY,FtL,{ier.)g^Df{e-ier.) (12) 

where G^^ and D^^ are the fermionic and bosonic Mat- 
subara Green's functions which are easily expressed in 
terms of the spectral functions and B^^q respectively. 
The Matsubara sums in the second lines of Eqs. and 
El only contribute to the real part of the self energies. 
The population factor is given in terms of Bose 

(&) and Fermi (/) population functions as. 



fuj — e 



(13) 
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We solved these equations numerically using bare Green's 
functions Gf^,[0,0], F^^j,[0,0] for calculating the self en- 
ergies Sf^. and $f;j- We show later that feedback effects 
give no significant changes within our model. 

Although we solve the equations above numerically 
without further approximations, some general remarks 
are in order. The function p^, ^ , , as a function of lo is 
at zero temperature nonzero only between w = and 
Lo — e, and is equal to sign(e) in this range. Because 
the spin fluctuation spectrum is gapped by much more 
than the thermal energy in the superconducting state, 
we can put for all practical reasons bn^^^ — 0. That 
means that we can neglect thermally excited modes, and 
only allow for emission processes at the resonant mode 
energy. For any gapped spin fluctuation spectrum with 
gap O, the first terms in Eg. 1111 and 1121 are negligible in 
the range —17 < e < (apart from temperature smear- 
ing near the value ±f2). Thus, assuming that the spin 
fluctuation spectrum is gapped below the resonance en- 
ergy, at zero temperature scattering of electronic exci- 
tations is disallowed in the interval — f2res < e < ^res- 
This is an expression of the fact that at least an energy 
fires must be spent in order to emit one spin fluctua- 
tion mode. This is the case for optimally and overdoped 
cuprates. For strongly underdoped cuprates, scattering 
is disallowed only in the range —Eg < e < Eg, where 
Eg is the spin gap which is smaller than fires- Also, as 
an implication, the renormalization function, determined 
by the real part of the self energy, is given in the low en- 
ergy range by the second terms of Fas. 1111 and 1121 onlv. 
In the following, we first consider the simple form of the 
model, which uses only the mode part of the spin fluctu- 
ation spectrum. After having gained some insight about 
the features caused by the resonance mode, we study the 
extended model which includes the continuum part as 
well. 



TABLE IL Minimal parameter set used in the calculations. 



Am 


^res 






g^WQ 


35 meV 


39 meV 


-34 meV 


2a 


0.4 eV^ 



We fit ARPES data near optimal doping,^ giving g^wg — 
OAeV^. This implies that the coupling constant is equal 
to g = 0.65 eV. This is a reliable value as discussed in 
Ref. [b^. In Table IHI we present our minimal parameter 
set entering the model (we only include the parameter 
from the band structure tight binding fit, as the results 
are insensitive to reasonable variations of the other pa- 
rameters), together with the values we used for optimally 
doped compounds. 



1. Electron Scattering 

We first discuss phase space restrictions for electron 
scattering in the d-wave superconducting state, and how 
they relate to the issue of whether the small relative 
weight of the resonance part of the spin fluctuation spec- 
trum leads to comparable in magnitude effects as in 
strong coupling superconductors. Using bare Green's 
functions, the self energy at zero temperature can be 
written as 



-^g'^WqAi^_gd{e + ilres + Ek-q) 

-Y,9''wqA+^^5[e~nres-Ek-q) (16) 



ReEf, 



E 



g'^Wq 



Ek- 



TT (fires + Ek-qY " 



(17) 



B. Contribution from the spin fluctuation mode 



For a sharp bosonic mode the spectral function is given 



by. 



g^B^^q 2g^Wq {S{UJ - fires) - S{llI + fires)) (14) 



where Wq is the energy integrated weight of the spin fluc- 
tuation mode, which is assumed to be enhanced at the 
Q = (7r,7r) point. Using the correlation length ^s/i, we 
write it as 



1 



4^2 



sfl 



(cos' 



2 M 



„2 Sjl ' 



(15) 



We will show below that it is a good approximation to 
assume the mode as perfectly sharp in energy, as cor- 
rections due to the finite energy width of the mode are 
negligible. From neutron scattering data obtained on 
Bi2Sr2CaCu208+5, the energy integrated weight of the 
resonance mode was determined as 1.9 /^^r^ leading (af- 
ter dividing out the matrix element 2/i^) to wq = 0.95. 



where Ek = y/^TJK^, = {l±^k/Ek)/2. The sum 
over q extends over the first Brillouin zone for the spin 
fluctuation momentum. For negative energies, only the 
first sum in Eq. 1161 is nonzero. The sum is a weighted 
average of the expression A'j^_^S{e + fires + Ek-q) with 
weight factors Wq. For given fermion energies, e, and 
momenta, fc, the delta function restricts the allowed spin 
fluctuation momenta q. Similar zero temperature formu- 
las hold for the off diagonal self energy. 



Im$^fc = -^g^WqCk-q 5{e - fires - Ek-q) 



-(5(e + fires + Ek-q) 



Re<i>5 ^ - y: 



g'^wq 

TT ' (f^res + Ek-q)'^ ~ 6^ 



(18) 
(19) 



with Ck = Ak/2Ek. 

In Fig. 13 we plot for k = kj^j and for several energies 
these restricted regions in g-space. The corresponding 
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FIG. 3: The black regions denote the part of the first Brillouin zone of the spin fluctuation momentum q which participates 
in scattering of electrons with momentum ku and energy e, as indicated above each picture. The amount of scattering events 
is controlled by the form factor for the resonance mode Wq, which takes its maximal value (wq ) in the center of the Brillouin 
zone at Q = Inside the black circles Wq > 'wq/2, and inside the white region Wq > wq/10. For small energies, a), 

only nodal electrons are scattered. For energies equal to Qres + Em, d), a large region around the A/-point of the fermionic 
zone participates in scattering events. Scattering electrons with this energy and momentum involves spin fluctuations with 
maximal weight, and thus almost exhausts the entire weight of the mode part of the spin fluctuation spectrum. Pictures c) 
and d) correspond to the special energies fires + A a and fires + Em, leading to cusp features in the energy dependence of the 
imaginary part of the self energy. 



weights for these regions, given by Wq, are maximal at (f = 
Q {lx = ly = tt), and decay away from that momentum. 
For reference, we define the regions inside the black circle, 
where Wq > wq/2, and the white regions, where Wq > 
wq/1Q. The calculations were done for finite T = 40K, 
and with a broadening parameter (5 = 5 meV in Eq. |31 

For energies —^res{= — 39meV) < e < 0, there is no 
phase space available for scattering. Scattering of elec- 
trons by the spin fluctuation mode sets in for e = — fi^es 
at q corresponding to the wavevectors q — {kM ~ kjq) 
mod (G), connecting the M point to the nodes (G de- 
notes a reciprocal lattice vector). In picture a) of Fig.O 
we show for e = —(fires + Jq^a) the mode wavevec- 
tors involved in scattering events. The weight for such 
events is very small, as can be seen from the fact that 
these wavevectors are outside the white region. Going 
further away from the chemical potential with e, the al- 
lowed mode wavevector regions increase, as shown in pic- 
ture b) for e = —{Vires + \^a)- When the special point 
e — —{Vires + A^) is reached (= —71.2 meV in our case), 
the arcs of g-regions involved in scattering events close at 
the points q — {ku — kA) mod (G), as shown in picture c), 
and electrons are scattered strongly between the M point 
and the A points. This leads to a cusp (or peak for very 
small quasiparticle broadening) in the energy dependence 
of the imaginary part of the self energy at this energy. 
Going further in energy, another special point is reached 
at e = -{VLres + En) (with Em = a/^^j + A^^), at which 
scattering events between the M points involving spin 
fluctuations with momentum q = Q (and with q = Q) are 
allowed. We show the corresponding regions in tf-space in 
picture d). This picture is important for understanding 
the large effect we obtain. First, the weight factor Wq is 



large in the patches of phase space for allowed scattering 
events around Q. Furthermore, because of the van Hove 
singularity in the band dispersion, these patches have a 
large area, almost filling the area inside the black circles 
in Fig. 13 This has as consequence that a large part of the 
weight of the resonance is exhausted for scattering elec- 
trons with energies equal to e = —{Vires + Em), which 
amounts to —87.8 meV for our parameter set. Going even 
further in energy, as shown in picture e), the amount of 
scattering events quickly decreases. The area which is 
involved in electron scattering events is maximal for en- 
ergies between 70 meV and 90 meV. For these energies, 
the involved spin fluctuations are also near the g-region 
where almost all their weight is concentrated. Thus, the 
strongest renormalization effects will take place in the 
energy range 70-90 meV. 

Let us compare this discussion with the case for con- 
ventional isotropic electron-phonon coupling. In this 
case, the weight factors Wq are constant. The relative 
amount of phonon wavevectors involved in scattering 
events is then equal to the ratio between the black areas 
shown in Fig. O and the total area of the Brillouin zone. 
This ratio is for the maximal case, picture d), equal to 
5%. That means that only 5% of the total phonon weight 
contributes to the imaginary part of the self energy. It is 
well known that electron phonon coupling easily leads to 
renormalization factors of the order of 2. In our case, the 
spin fluctuation weight of the mode is only about 5% of 
the total spin fluctuation weight, but it is concentrated 
in the region inside the black circles in Fig. O Almost 
the total area inside the black circle contributes in the 
case of picture d) , showing that the same amount of only 
a few percent of the bosonic spectrum is involved as well 



8 



e-(^^res+EM) 












FIG. 4: The same as Fig. El d) for a fermionic wavevector 
at the nodal point, k — kN- Because the aUowed region for 
scattering events is outside the region of enhanced spin fluctu- 
ations, the corresponding cusp feature in the imaginary part 
of the self energy is weaker than for electrons with momenta 
near the M-point. 



for spin fluctuations in high Tc cuprates as for phonons in 
conventional strong coupling superconductors. Thus, the 
renormalization of the fermionic dispersion is expected to 
be of the same order of magnitude, and our explicit cal- 
culations confirm this. 

In Fig. ^ we show the (f-space areas corresponding to 
Fig. Old), but for electrons near the nodal wavevector. As 
can be seen, the feature due to the van Hove singularity 
region is now weighted by a smaller value of Wq . Because 
of this, for nodal electrons, the corresponding peak in the 
self energy is smaller than for momenta near the M point. 
It turns out that for the nodal electrons, the feature at 
— {Qres + A^) is more pronounced than that at —(fires + 



Em). 



2. Renormalization factor and electron lifetime 

The self energy has a characteristic shape as a func- 
tion of energy, which is conserved qualitatively for all 
points in the Brillouin zone. This is a consequence of 
the fact that all points are coupled via the spin fluctu- 
ation mode, which has a finite width in momentum, to 
all special points in the Brillouin zone with their cor- 
responding characteristic energies. These special points 
are the nodal N points, and the van Hove singularities 
at the M points and the A points (the latter is a disper- 
sion maximum in the superconducting state). Because 
the general shape of the energy dependence of the self 
energy does not vary much with momentum (although 
the overall intensity does), it is sufficient to discuss the 
important features in the energy dependence of the self 
energy at the M point. 

We numerically evaluated the self energy, using a 
broadening parameter 5 = 5 meV. In Fig. |S1 we show 
the results for the renormalization function and electron 
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FIG. 5: Renormalization factor at the M point, Zmie) (top) 
and electron scattering rate at the M point, rAf(e) (bottom). 
The thin lines denote some characteristic energies: if^res 
(dotted), ± (fires + Aa) (dashed), and ±(r2res + Em) (dot- 
dashed). Electrons at low temperatures are scattered only 
if their energy is larger than Qres, so that they are able to 
emit a collective mode excitation. The parameters used are: 
T = 40K, fires =39 meV, Am =35 meV. 



scattering rate at the M point, 

ZM{e) = l- ^""^^'^^'^ Tm (e) = -Im SJ:,, (e) (20) 
e 

as a function of energy. 

There are three characteristic energies (in addition to 
temperature, which smears all features by fc^T). Region 
I is bounded by the resonance energy, fires, and has zero 
scattering rate at zero temperature (this statement is true 
for electrons at any point in the Brillouin zone). At fi- 
nite temperature, a region fcsT around iiflres allows for 
a small amount of scattering, even in region I. Because 
states are occupied near the M point, we will only dis- 
cuss negative energies in the following. At e = —fires, 
scattering for all electrons in the Brillouin zone sets in 
due to coupling to nodal electrons via emission of a spin 
fluctuation mode. Absorption processes are negligible 
due to the large (compared to temperature) mode en- 
ergy. In region II, a larger and larger area around the 
nodes participates in scattering events, (as can be seen 
from pictures a and b in Fig.O, until finally the point at 
the zone boundary with maximal gap, ±A^, is reached 
(picture c in Fig.jSJ. This point corresponds in Fig.|31to 
a cusp feature in the imaginary part of the self energy at 
-(flres+^A)- The third feature, at —{flres+EM), corre- 
sponds to the van Hove singularity at the M point of the 
Brillouin zone, which is close to the chemical potential 
in cuprates (picture d in Fig. O . The proximity of this 
van Hove singularity leads to a stronger peaked feature 
in the scattering rate near ±(flres + ^a) compared to the 
case where this van Hove singularity at the M point is 
absent. The renormalization factor is rather constant in 
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region I as a consequence of its connection to the imag- 
inary part via Kramers-Kronig relations. The enhance- 
ment in regions I and II compared to unity comes from 
two step features at ±(f2res + Em) and at ±(f2res + A^). 
Note that the step feature due to the van Hove singular- 
ity at the M point contributes about 50% to the total 
enhancement. The small features at ±i7res are due to 
the finite lifetime of the electrons involved in scattering 
processes as discussed below. The onset of scattering at 
the emission edge for the spin fluctuation mode occurs 
as a jump if the electrons involved have a finite spectral 
width. At even higher energies, corresponding to Fig. |2| 
e), the scattering due to the spin fluctuation mode be- 
comes less effective. Note that the spectral peak of the 
electrons at ±Afc is either in region I or in region II. Thus, 
quasiparticles near the nodal regions are always sharper 
in energy then quasiparticles near the maximal gap re- 
gions. In overdopcd cuprates, the maximal gap is usually 
smaller than the mode energy, so that for the broadening 
of the quasiparticle peaks, the spin fluctuation mode is 
not relevant. 

For the following discussion, it is useful to derive ap- 
proximate analytical expressions. At zero temperature, 
using Eq. 1171 we obtain 



^Af(e) = 1 + 



1 



9 



{fir 



Ek^ 



(21) 



The main contribution comes from the regions where 
Ek^f-q is less than 100 meV. We can estimate those re- 
gions by the requirement that fcjv/— (fis in the area around 
the M points deliminated by ±0.357r in M — 1^ direction 
and by about O.Stt along the M — T direction. Then, 
replacing I^kM-q by "^m, and Ek^^-q by Em, we per- 
form the (f-sum over that area of the function Wq. We 
denote X^n^? over this area by Iq. For our model we 



have /q = 
Zm{<^) 

AN) 



0.035. Using this approximation, we obtain 
1 



1 



EmY 



^ + Aif(e) (22) 



Here, \\j (e) denotes the contributions coming from the 
regions where k' = k^i — q outside of the above range. 
It is dominated by contributions where k' is near the 
nodal regions of the Brillouin zone, thus the relevant spin 
fluctuation momentum is q ^ (/cm — k^) mod G. The 
contribution Aj^^'' is smaller than the first term in Eg. 1221 
but not negligible. Because Eq. 122 neglects the disper- 
sion between (fires + A^) and {fires + Em) near the M 
point, it should be used for energies not too close to the 
region between these two values. We will make use of 
this formula below for energies near e = Em, where this 
formula gives a good approximation. 



3. The quasiparticle scattering rate 

For overdopcd materials the quasiparticle peak at the 
AI point is situated below the onset of scattering due to 
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FIG. 6: Renormalization factor at the M point, ZM{e) (top) 
and electron scattering rate at the M point, rAf(e) (bottom). 
The picture compares results for two different residual quasi- 
particle linewidths: 5 = 1 meV (full lines) and 5 = 5 meV 
(dashed lines). As dotted lines the nodal contributions, when 
restricting the quasiparticle momenta to the regions outside 
the area around the M points discussed in the text, is shown 
for (5=1 meV. 



emission of spin fluctuations. In this case the width is de- 
termined by other processes, and we model this residual 
quasiparticle width by a parameter b. In Fig. |3 we show 
the influence on the renormalization factor and the scat- 
tering function of the residual quasiparticle width. We 
compare the results for (5 = 5 meV with those for (5=1 
meV. For very small quasiparticle broadening (full lines) 
the cusp features in the imaginary part of the self energy 
turn into peaks (which ultimately evolve into square root 
singularities for perfectly sharp quasiparticles and reso- 
nance). The second feature to mention is that the scat- 
tering rate near the onset points, ±f2res, is influenced 
strongly by the residual quasiparticle width. Because 
this onset region governs the quasiparticle width in un- 
derdoped cuprates, as we show later, we study it in the 
following in more detail. In the lower part of Fig. |B1 we 
show as a dotted line the contribution to the electron 
scattering rate coming from the final states not too close 
to the M points (the regions which determine A^^\ in- 
troduced above) as compared to the full scattering rate 
(full line) . It is clearly seen that the sharp features come 
from the M point regions, whereas the nodal regions con- 
tribute to the onset of electron scattering and provide a 
smooth constant background at higher energies. 

The behavior of the imaginary part of the self energy 
near the onset points, ifires, in Figs. Eland El is deter- 
mined by the nodal electrons. For larger residual quasi- 
particle widths ((5 = 5 meV, dashed lines in Fig. jSJ there 
are states available at the chemical potential (coming e.g. 
from impurity scattering), which increase the number of 
final states for scattering events. Thus, the onset in Fig. 
El for the electron scattering rate is stronger in this case 
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than for 5=1 meV. For zero temperature there wiU be 
a jump at energy iVlres in the imaginary part of the self 
energy, which causes the small cusps at the same energy 
in the renormalization factor (top panel in Fig. 1^. For 
(5 = the onset is linear in energy. 

We will estimate analytically the onset behavior near 
these points for the case 6 = now. For this we use 
Eas. 1101 and [TCI We replace qhy k' — k, approximate Wq 
by wmn — WkM-kNJ s-nd linearize the dispersion around 
the nodes, = v^ik' — fcAr), (,k = VNik' — kpf). Here 
VA = dk^k and vn = dk£,k taken at the N point. For our 
model, we have t^A = Am sin(fca;Af)/\/2, which is valid 
near optimal doping. (But note that for underdoped 
cuprates, wa was experimentally shown to be smaller 
than that value, perhaps scaling with ksT^ instead of 
with Ajvf •^*) Performing the k' sum and summing over 
all four nodes, we arrive at 



Im<5EW(6) 



g WMN 



\e + r2res|B(-e - fires) 



+ \€ - nres\Q[i - ^res)) (23) 

Here, the 0-f unction is unity for positive argument and 
zero otherwise. Thus, the slope of the scattering rate 
at e = ±Ores is given by T9'^wmn /t^vnva- For the 
parameters in Tables and m the magnitude of this 
slope is equal to 9.5 wmn/wq « 0.56. Note that Eq. [53 
gives a good approximation of the scattering rate in the 



interval Vlr 



< lei < rir 



Aa/2. For energies fur- 



ther away from the onset, the change of the quantity va 
(which goes to zero at the A point) leads to a stronger 
increase. Finally, for underdoped cuprates the excita- 
tion energy at the M point, Em, is larger than fires- 
Then, the quasiparticle linewidth at the M point is given 

by —lm.6T,^^\—EM)/ZM{—EM)- Thus, for underdoped 
cuprates it is given by. 



M 



g'^WMN Em — 



Z 



M 



(24) 



with Zm = Zm{~Em)- Near the nodes, on the con- 
trary, the quasiparticles will stay relatively sharp even in 
underdoped compounds because the peaks positions are 
then below the onset energy ±r2res- 



4- The coupling constant and the weight of the spin 
resonance 

One potential criticism of a model which assigns the 
observed anomalies in the dispersion to coupling of elec- 
trons to the spin resonance mode is the spectral weight 
of the resonance, /q, which amounts to only a few per- 
cent of the local moment sum ruleiS^ Our calculations 
show that this is not an obstacle^^ as we obtain a di- 
mensionless coupling constant of order one, as observed 
experimentally. 



Here we estimate \m, given by Zm{^) — 1, for the 
resonance mode. From Eq. it is equal to 



A 



M 



1 



(f7r 



E 



M 



- + Aif(0) (25) 



Using values for optimal doping f Table Hll). the first term 
in this sum is equal to 17.44/o, which amounts to about 
0.61 (in our model /q = 0.035). This is already a large 
part of the total coupling constant, which from Fig. [S] 

is \m ~ 0.9. The contribution 1 + X^^\e) is shown as 

dotted line in the upper part of Fig. El A^^"* is not neg- 
ligible, but contributes about 30% to the total coupling 
constant. 

We obtain an analytic formula for the low energy cor- 
rection to the renormalization factor due to scattering be- 
tween nodal points and M points, X^j^\e), by a Kramers- 
Kronig transform of ImSj^^^ (e) , in which only energies up 
to a cut-off zL{flres + A^) are taken into account, and re- 
placing ImEjj^''(e) above this cut-off by a constant (see 
the dotted lines in Fig. |SJ) equal to its value at the cut-off. 
The result for e = is, 



Aif(O) 



g WMN 2 



-In 1 



A, 



(26) 



For our parameter set this amounts to A^^-'(O) ~ 0.21. 

Note that A^^ increases with decreasing fires- 

To summarize, dimensionless coupling constants (com- 
parable to those for strong-coupling electron-phonon sys- 
tems) are easily achieved with reasonable parameters by 
coupling electrons to the spin resonance. 



5. Particle hole asymmetric renormalizations 

From Eq.EI we see that the second term in the numer- 
ator, proportional to ^k-qi affects the band dispersion ^k- 
The resulting renormalization is given by. 



£,f.,k = Cfe 



9 Wq 



k — q 



{fires + Ek-q)^ - 



(27) 



From this formula, it is clear that notable renormaliza- 
tions of the Fermi surface only take place if ^k-q is not too 
far from (but also not at) the chemical potential. Thus, 
the largest renormalizations are expected at the M point 
regions of the Brillouin zone. 

In Fig. [7| (left), we show the particle- hole asymmetric 
part of the self energy as a function of e for electrons at 
the M point of the Brillouin zone. The imaginary part 
shows a peak due to the van Hove singularity at the M 
point, but the cusp feature due to the A points is missing, 
because points where k — q is on the Fermi surface do not 
contribute to the sum in Eq. 1271 The real part indicates 
that the renormalization of the dispersion is confined to 
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FIG. 7: Left: The part of the self energy (eV) defining the 
band renormahzation at the M point as a function of en- 
ergy, for Am=35 meV, ^lre3~39 meV, and g=0.65 eV. Right: 
The quantity ^A/(e)/(^Af (e)^*/) (filled circles) and 1/Zm{^) 
(empty circles) for e — —Am as a function of the coupling 
constant g. 



energies between —ilres — Em and fires + Em- Using the 
same approximation procedure as above, we obtain for 
the renormahzation at the M point. 



ilres + Em 



Em i^re 



EmY 



(28) 



The first important point is that the renormahzation has 
opposite sign to ^fc, thus the band is renormahzed to- 
wards the chemical potential. In particular, there is a 
'pinning' effect of the van Hove singularity at the M point 
to the chemical potential, as long as is of the order 
of r^res- Furthermore, the renormahzation factor Zm{^) 
from Eq. 122 increases this effect, as £,m/Zm defines the 
quasiparticle dispersion. 

In order to quantify this, we show in the right 
panel of Fig. [7| the relative changes of the dispersion, 
CM(e)/(^M(e)CA/) (filled circles), in comparison to the 
inverse renormahzation factor l/ZM{f) (empty circles). 
The latter would give the band renormahzation in the 
absence of particle hole asymmetric parts in the self en- 
ergy. As can be seen in this figure, the band is renormal- 
ized towards the chemical potential and even crosses it 
for large coupling constants. For coupling constants near 
0.6 eV, the renormahzed band is close to the chemical 
potential. Thus, the dispersion of the peak in ARPES is 
negligible in the M point regions as a result of the renor- 
mahzation of the dispersion. The renormahzation of the 
band implies an increase in the chemical potential, so as 
to keep the particle density constant. This effect would 
increase the distance between the chemical potential and 
the van Hove singularity at the M point, leading to an 
equilibrium value in a self consistency loop. We did not 
solve this self consistency problem, but assumed that our 
parameter choice is close enough to the self consistent so- 
lution to capture the main physics. 



6. Off diagonal self energy 

In order to understand the renormahzation of the order 
parameter due to coupling to the resonance mode, we 



observe from Eq. ^| 



A,,fc = Afe 



(29) 



This formula is very similar to that for the band renor- 
mahzation, except that the order parameter at momen- 
tum k — q now determines the renormahzation effect. 
Note that if Wq were independent of q, no renormahza- 
tion would take place due to the d-wave symmetry of the 
order parameter. Since the spin fluctuation continuum, 
which we discuss later, is very broad in momentum, the 
renormahzation effects in the off-diagonal components is 
dominated by the resonance contribution. As the order 
parameter vanishes at the node, we concentrate on the 
renormahzation near the M point region again. Adopt- 
ing the approximations as above (note that contributions 
from the nodal regions cancel because of the d-wave sym- 
metry) , and using Eq. 1191 we arrive at 



AM(e) 



1 



E 



M 



TT Em {^res + Em)^ ~ 



(30) 



The positive sign is due to the fact that Am+q — — Ajv/. 
As a result of this, there will be a compensating effect 
when calculating the quantity AM(e)/-^Af(e), which de- 
termines the peak position. In Fig. |H1 (left), the real 
and imaginary parts of the off diagonal self energy at 
the M point are shown. The imaginary part is relevant 
only for energies with absolute value > ilres + A^. For 
smaller energies, the main effect is to increase the mag- 
nitude of the order parameter A^ in the energy range 
—fires — Ayi < e < fires + ^A- Notc that the self en- 
ergy due to coupling to the resonance mode has d-wave 
symmetry, like the order parameter. Thus, the coupling 
to the resonance mode supports superconductivity. In 
order to quantify the amount that the resonance mode 
contributes to the spectral gap, we show in Figs. |S1 (right) 
and the quantity AM(e)/(^A/(e)AM) (together with 
1/Zjvf (e) for comparison) as a function of three different 
parameters: g, Am, and fires- 

As can be seen from these figures, although the renor- 
mahzation factor Zm would reduce the order parameter 
considerably, the off diagonal contribution to the self en- 
ergy from coupling of electrons to the resonance mode 
restores the gap to its original value. Thus, the reso- 
nance contribution to the gap is as big as that from other 
sources, and starts to dominate if the coupling constant 
exceeds about 0.5 eV. 

The reason why Am(Am) is so close to Zm{Am)Am 
is that the additional factor 1 -I- fires/ Em in Eg . 1301 com- 
pared to Ea. l22l is approximately canceled by the presence 
of the additional A^^'' (Am) in Eq.|221 An analogous term 
in Eq. 1201 is missing due to the sign change of the order 
parameter at the node. The degree to which this cancel- 
lation holds is a surprising numerical result and allows us 
to avoid a self consistency loop for the determination of 
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FIG. 8: Left: The off diagonal self energy (eV) at the 
M point of the Brillouin zone as a function of e is shown 
for coupling constant (7=0.65 eV. Right; The quantities 
AA/(e)/(^M(e)AM) (filled circles) and 1/Zm{^) (empty cir- 
cles) are shown for e = — A^f as a function of the coupling 
constant g. Parameters used are Am=35 meV, Slres=39 meV. 
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FIG. 9: The quantities AM(e)/(^Af(e)AAf) (filled circles) 
and l/ZM(e) (empty circles) are shown for e = —Am as a 
function of Am (left, for g = 0.65 eV and Qres = 39 meV) 
and fires (right, for g = 0.65 eV and Am = 35 meV). 

Am near optimal doping. Thus, the experimental param- 
eters which enter our calculations are already sufRciently 
self consistent. 

7. Spectral functions at the M point 

In this part, we discuss the spectral lineshape, which is 
an experimentally accessible quantity. The main features 
of the spectral lineshape are captured in the simple model 
neglecting the continuum part of the bosonic spectrum. 
We discuss in the following the influence of the different 
parameters of the theory on the spectral function, 

A(e,fcM) = -2ImG«(e,fcM) (31) 

and will discuss changes due to the continuum part of 
the spin fluctuation spectrum later. In our numerical 
studies, we used a broadening parameter S — 5 meV. 
This accounts for processes not covered by scattering by 
spin fluctuations. 

In Fig. EH we present the results for the spectral func- 
tion at the M point of the Brillouin zone for both a per- 
fectly sharp resonance and for a finite width of the reso- 
nance of 10 meV. It is obvious that the energy width of 
the resonance has very little effect on the ARPES spec- 
tra, except a slight reduction of the peak height. 
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FIG. 10: Spectral functions at M for a perfectly sharp reso- 
nance (full line) and for a resonance with a finite energy width 
of 10 meV (dashed line). Parameters are for optimal doping. 
The finite width of the mode has very little influence on the 
ARPES spectra, and can be neglected for most purposes. 



Thus, we will concentrate all our following discussions 
on a perfectly sharp resonance mode. The main features 
of the spectral function is the dip feature at an energy 
of about the resonance energy relative to the peak. The 
peak position at ~Em is renormalized by self energy ef- 
fects discussed above, and is shifted from the bare —Em 
to be near —Am- The dip feature is actually spread out 
over a range of size Em — A a, and it is the onset of this 
dip feature which defines the resonance energy, fires ■ The 
dip feature is followed by a hump at higher binding en- 
ergies, and the position of the hump maximum is very 
sensitive to the coupling constant and to damping due to 
the spin fluctuation continuum, as we show later. Thus, 
we concentrate in the following on the peak-dip structure. 
Another feature worth mentioning is the asymmetry of 
the lineshape at positive and negative binding energies, 
with a relatively weak dip feature on the unoccupied side 
compared to the occupied side. 

In Fig. ^2 (left), the effect of a varying resonance en- 
ergy Qres (keeping all other parameters at their values for 
optimal doping) is shown. The spectral function shows 
two effects. First, the peak weight is reduced with de- 
creasing mode energy. Second, as soon as the quasiparti- 
cle excitation energy exceeds fires, strong damping sets 
in. We can understand these results in the light of the 
discussion for the self energy. As we mentioned above, 
the scattering rate has a gap equal to fires- Thus, as 
long as the spectral peak is situated below that energy, 
in region I of Figs. 1101 and there will be no damping, 
and the peak width is set by the residual broadening due 
to other processes. If the peak is positioned above fires 
(region II in Fig. 110(1 . it feels the self energy in region II 
of Fig. 13 and will be broadened. Because in region II the 
self energy is dominated by scattering processes involv- 
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FIG. 11: Spectral functions at the M point for varying fires 
(left, for 10 meV, 20 meV, 30 meV and 40 meV from bottom 
to top; the thin lines denote the value ±Am), and for varying 
Am (right, for 15 meV, 25 meV, 35 meV, 45 meV, 55 meV, 
and 65 meV from bottom to top; the thin lines denote the 
value ±Qres)- All other parameters are kept fixed at their 
optimal doping values. The spectra are offset for clarity. 



However, in this case, the weight of the peak is affected 
only weakly, as we will discuss below. 

Finally, we show in Fig. 1121 the influence of increas- 
ing coupling g, and of an increasing distance of the van 
Hove singularity from the chemical potential, ^a/- In 
both cases, the hump energy is strongly affected, moving 
to higher binding energy with increasing coupling and 
increasing ^m- In the left panel, one can also see that 
the weight of the peak is strongly reduced with increas- 
ing coupling constant. This is not the case with varying 
£,M , as seen from the right panel in Fig. ^1 and will be 
discussed in more detail below. 



The coherent quasiparticle weight of the ARPES 
spectrum 
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FIG. 12: Spectral functions at the M point for varying 
coupling constant (left, g = 0.1 eV, 0.3 eV, 0.5 eV, 0.7 eV, 
and 0.9 eV from top to bottom; thin lines denote ±Am) and 
for varying distance of the van Hove singularity at the M 
point from the chemical potential (right, for — ^aj = meV, 
20 meV, 40 meV, 60 meV, and 100 meV from top to bottom; 
thin lines denote ±AAf). All other parameters are kept fixed 
at their optimal doping values. The spectra are offset for 
clarity. 

ing nodal electrons, the width in this region is set by the 
imaginary part of the self energy divided by the renor- 
malization factor, and is given in Eq. |^ At the same 
time, for decreasing resonance mode energy, the incoher- 
ent part of the spectral function grows, taking weight 
from the quasiparticle peak. 

Thus, in Fig. ^2 which is for Aat =35 meV, the 
quasiparticle weight increases from the lowest curve (for 
fires = 10 meV) to the uppermost curve (for fJ^es = 
40 meV). Simultaneously the broadening decreases. As 
the onset of quasiparticle damping and the loss of the 
coherent part of the spectrum is a result of a decreas- 
ing resonance mode energy relative to the gap, the same 
effect is expected by increasing the gap keeping the res- 
onance mode energy constant. This is shown in Fig. 1111 
(right). In this case, the onset of quasiparticle damp- 
ing is always at the same energy fi^es =39 meV, but 
for the lowest curve, corresponding to a small gap of 15 
meV, quasiparticle peaks are well established, whereas 
for the uppermost curve, corresponding to a large gap of 
65 meV, the quasiparticle peaks are strongly broadened. 



Although one can define a quasiparticle residue via 
the renormalization factor Z{e), in light of the experi- 
mental studies, we will in this part study the weight of 
the quasiparticle peak in the ARPES spectrum, deter- 
mined by numerically integrating over the peak region. 
For strongly renormalized spectra, this experimentally 
motivated quantity will differ from the first. We note 
that due to coupling to the mode, the peak weight is re- 
duced and redistributed to the hump. Because the peak 
weight in the experimental literature is often referred to 
as the 'coherent quasiparticle weight', we will use the 
same terminology here. 

We consider the spectral function at the M point of 
the Brillouin zone. Because the peak is separated from 
the hump by a dip which extends from — ei = — (rj^es + 
I^a) to — e2 = — (f^res + Em), we define as the coherent 
quasiparticle weight the quantity, 

ZM = -- I deImG^,(e) (32) 

Without interactions between the quasiparticles, z — 0.5 
at the Fermi surface, because the quasiparticle peaks at 
±A in BCS theory each have one half of the total weight; 
the value at negative energy is somewhat larger than 0.5 
at the M point because it is an occupied state. Coupling 
of the quasiparticles to the mode reduces z. In Figs. [T^ 
and 1141 our numerical studies are summarized. The re- 
sults are as follows: (1) zm is only weakly dependent on 
the gap and the band structure in the relevant parameter 
range; (2) za/ is proportional to the mode energy fires! 
together with the experimental finding fi^es oc ksTc, this 
means za/ oc ksTc] (3) for coupling constants of order 
the band width or larger, zm c>c 1/ [g^wq); for smaller 
coupling constants, 1/ zm ^ A + Bg^wq with A and B 
constants; (4) zm weakly decreases with increasing anti- 
ferromagnetic correlation length ^^fi . We can understand 
some of these features using the approximate expression 
of Eq. 123 Evaluating Zm{£) at e = —Em, and tak- 
ing into account the coherence factor at the M point, 
Ajj = A^j(—Em), and the nodal renormalization factor 
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FIG. 13: The coherent quasiparticle weight as a function of 
Am (left) and (right) for Qres =39 meV. Although the 
peak width changes considerably as a function of Am, the 
peak weight is only weakly dependent on Am and ^m- 




FIG. 14: The inverse of the coherent quasiparticle weight 
1/z is approximately a linear function of g^wq (left). Here 
we have chosen Q,res =39 meV and Am =35 meV. The right 
panel shows that z is proportional to fires- 



W^l + AW(-i?M), gives 



ZM 



O A~ 



(33) 



which defines the constants A and B. 

In the underdoped region, where ^res is much smaller 
than 2Em, we can approximate further to obtain 



Wo 



(34) 



Here, we neglected the first term in the denominator of 
Eq. 1331 compared to the second, which is justified when 
zm is small. In the overdoped region, where g^Io de- 
creases and ^Ires approaches 2Ah (where A^i is the gap 
at the hot spots), this scahng with fi^es should break 
down according to Eq. Note that experimentally, the 
relation fi^es ~ -i-^ksTc was shown, and also the rela- 
tion 



zmEm \ 

/ (exp) 



ksTc 



0.5 



(35) 



was experimentally found»^ Thus, our expression Eg. 1341 
would be consistent with the experimental finding if with 
doping E\j scaled with g^Io- Within our theory this ex- 
perimental finding can be interpreted as an indication 



that the phenomenological order parameter is gov- 
erned by the same coupling constant g. 



C. Contribution of the spin fluctuation continuum 

At energies higher than that corresponding to the con- 
tinuum edge of the spin fiuctuation spectrum, additional 
broadening due to coupling to that part of the spectrum 
sets in. Because the continuum extends to electronic en- 
ergies (~ eV), the introduced scattering rate will increase 
continuously with energy up to electronic energies as well. 
We model the continuum part by 

g^Bl^^ = 2g^c^ - 2A,0 - e(-c^ - 2A0) (36) 

where the gap in the continuum spectrum is given by 
2A^. This form for the gapped continuum is similar to 
the gapped marginal Fermi liquid spectrum considered 
earlier by other authors .-^s^S. The momentum dependence 
takes into account the experimentally observed flatter be- 
havior around the (tt, tt) wavevector at higher energies, 
and is modeled as 



1 + (32C,4)- 



1 



16^4 (cos4 ^ 



(324')- 



C0s4f ) 

(37) 

with a correlation length ^ = 0.5a compatible with ex- 
perimental findings. We subtracted a background term, 
so that the response far away from the (tt, tt) wavevector 
is small, as experimentally observed (we have chosen this 
background term so that Cq is zero at q = 0). 

For the chosen correlation length, the momentum av- 
erage of Cq gives 0.5cQ. The constant cq can be obtained 
from the experimental values for the momentum aver- 
aged susceptibility at 65 meV, which was found to be 
6/x^/eV for underdoped YBa2Cu307_5 in the odd chan- 
nel, and about 3/i^/eV in the even channelj^ Dividing 
out the matrix element 2/i^, this gives cq « 6/eV and 
3/eV respectively. The corresponding values near opti- 
mal doping should be smaller. We use in our calcula- 
tions Cq — 5.6/eV and g = 0.65 eV. The choice of this 
value is motivated by the ARPES measurements on opti- 
mally doped Bi2Sr2CaCu208-5 of the high energy (hnear 
in excitation energy) part of the momentum linewidth, 
which gives = 0.75e i^Tii? This coupHng includes both 
the even and odd (with respect to the bilayer indices) 
contributions of the spin fluctuations, in contrast to the 
coupling to the mode, which is present only in the odd 
channel. Note that our value for cq is about a factor 1.6 
smaller than neutron scattering measurements give for 
underdoped YBa2Cu307_5. Because in optimally doped 
compounds the intensity of the spin fluctuation contin- 
uum is smaller than in underdoped ones, this is a reason- 
able value for optimal doped Bi2Sr2CaCu208-5 

The spin fluctuation continuum is gapped in the odd 
channel from zero energy to twice the gap at the 'hot 
spots', 2A/i, which is slightly less than twice the maximal 
gap. This means that additional damping only sets in for 
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FIG. 15: The different contributions to the renormalization 
factor (top) and the scattering rate (bottom) are shown for 
the M point (left) and for the A'^ (node) point (right). Dotted 
curves are the contribution from the spin fluctuation contin- 
uum, dashed the contribution from the spin fluctuation mode, 
and fuU both contributions. 

energies |e| > 2A;i. This corresponds in optimally doped 
compounds to about 65 meV. In the even channel, the 
optical gap (~60 meV) persists into the normal statei^ 
The continuum formally has to be cut-ofF at high en- 
ergies. This cut-off does not affect the imaginary part of 
the self energy, but its choice leaves a real term of the 
form — Ce at energies small compared to the cut-off en- 
ergy scale. This term, equivalent to a contribution to the 
renormalization factor which is constant up to the high 
energies, has to be regarded as an additional phenomeno- 
logical parameter. The constant C depends on the model 
one uses for the high energy tail of the spin fluctuation 
spectrum. Because we model the continuum by a con- 
stant, which overweights high energies, we have chosen a 
relatively low cut-off of 200 meV for our model spectrum. 
Because the constant C is only weakly (logarithmically) 
dependent on the cut-off, the exact energy of the cut-off 
is not crucial. 

In the simple form of our model, we absorbed the 
renormalization from the continuum into the band dis- 
persion ffc. Now we take into account explicitly the con- 
tinuum, and thus have to start with a band dispersion not 
renormalized by this contribution. We found that we can 
reproduce experiment best by rescaling the dispersion 
from Table^in the following way: ^^"'^'"^ = 1.5^fc— 0.5^jv/- 
With this choice, the van Hove singularity at the M point 
has the same distance from the chemical potential as be- 
fore. 

In Fig. El the continuum contribution to the self- 
energy is shown as a dotted line. As can be seen from 
the figure, the continuum contribution to the scattering 
rate sets in above the structures which are induced by the 
mode. It also contributes considerably to the renormal- 
ization factor. As mentioned above, the renormalization 
does not decay up to energies of 200 meV, consistent with 
experiment. At the nodal point, the modification due to 
the continuum relative to the mode part is strongest. The 
importance of the continuum contribution can be seen by 
noting the strong similarity of the lower right hand panel 
of Fig. El to self-energies extracted from ARPES data 
along the nodal direction>iSi2S 



Finally, note that in the normal state, the even chan- 
nel stays gapped. That means that at the N point, the 
self energy for scattering between bonding bands and be- 
tween antibonding bands (but not between bonding and 
antibonding) is similar to one half the continuum contri- 
bution (dotted line) in the right panel of Fig. El This will 
induce a weaker kink feature in the normal state at an 
energy equal to the even channel (optical) gap in the spin 
susceptibility, which is around 50-60 meV. Correspond- 
ingly, the high energy renormalization will be present in 
the normal state, but weaker. The difference between 
the high energy renormalization in the normal and su- 
perconducting states is mainly due to the appearance of 
a continuum gap in the odd channel. The low energy 
renormalization is mainly due to the appearance of the 
mode in the odd channel. 



D. Renormalization of EDC and MDC dispersions 

In the following, we discuss the dispersion of the spec- 
tral lineshape through the Brillouin zone and study the 
corresponding EDC (as determined from the spectral 
maximum as a function of energy) and MDC (as de- 
termined from the spectral maximum as a function of 
momentum) dispersions. We include both the mode and 
the gapped continuum of the spin fiuctuation spectrum. 
In Figs. I16I19I we show dispersions of the ARPES spec- 
tra along several selected paths in the Brillouin zone. In 
the left panels of the figures, the intensities and spectral 
lineshapes can be followed, and in the right panels, the 
corresponding dispersions of the peak maxima and hump 
maxima in the EDCs are shown as circles, and the max- 
ima in the corresponding MDC dispersions as curves. A 
general remark concerns the linewidth of the high energy 
features compared to the low energy features. Due to 
the strong self energy damping effects setting in above 
the dip energy (Fig. I15|) , the hump features are consider- 
ably broader than the peak features for all momenta in 
the Brillouin zone. This holds for both EDC and MDC 
dispersions. Note that even without taking into account 
the lifetime effects due to the spin fluctuation continuum, 
the high energy features are much broader in energy than 
the low energy features^ To account for the experimental 
MDC linewidth, however, one has to take into account 
the continuum contribution. 

Starting with Fig. El we follow the dispersion along 
a cut going from the M point of the Brillouin zone to- 
wards the Y point. The A point corresponds to spectra 
roughly in the middle of the set. From the left panel, 
we see that sharp peaks are restricted to the momentum 
regions between the M and A points. The dip structure 
is maximal at the M point and much weaker at the A 
point. The corresponding dispersion, shown in the right 
panel, reproduces the experimental findings^^ of two al- 
most dispersionless EDC branches, one for the peak and 
one for the hump. The MDC follows the peak branch, 
then shows a nontrivial variation at energies within the 
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FIG. 16: Left: Dispersion of the spectral intensity and line- 
shape as a function of momentum along the M — Y cut, 
{ky = n, — 0...0.47r in steps of 0.047r from top to bottom). 
Right: EDC (circles) and MDC (curve) dispersions from max- 
ima of the curves shown in the left panel. In the EDC disper- 
sion, the low energy peak and the high energy hump with the 
break feature in between is clearly visible. Because the bot- 
tom of the normal state dispersion is at ~ — 34meV, the 
MDC shows only a broad maximum at Af for high energies. 




FIG. 17: Left: Dispersion of the spectral intensity and line- 
shape as a function of momentum along the M—F cut {kx = 0, 
ky = 0.67r..7r in steps of 0.047r from top to bottom). Right: 
EDC (circles) and MDC (curve) dispersions from maxima of 
the curves shown in the left panel. 

gap edge. This behavior is discussed in Ref. |U The 
Fermi crossing is only slightly shifted with respect to the 
unrenormahzed value of kx — O.IStt. At higher energies, 
the MDC is peaked at M. 

Going from the M point in the direction of the F point, 
the corresponding dispersion of the ARPES spectra is 
shown in Fig. El On the left side, one can see that 
the intensity of both the peak and the hump is almost 
unaffected in the region between the M point and roughly 
O.Stt from there in direction of F. In this range, the 
renormalized EDC dispersion of the hump is extremely 
flat, as seen in the right panel, and the peak shows a 
moderate dispersion, becoming almost flat between qy = 
0.97r and qy =^ tt. When going further away from the 
M point, the intensity of the peak drops sharply, and a 
strong dispersion of the hump sets in. There is a clear 
break between the peak and the hump EDC dispersion 
due to the dip. The MDC along this cut follows the 
peak near the M point, but changes over to the hump 
dispersion at roughly the point where the hump starts 
to disperse strongly away from the chemical potential. 



FIG. 18: Left: Dispersion of the spectral intensity and line- 
shape as a function of momentum ky = O.Bvr, k^ = 0...0.47r in 
steps of 0.027r from top to bottom. Right: EDC (circles) and 
MDC (curve) dispersions from maxima of the curves shown 
in the left panel. 



In this range, at energies between 70 meV and 100 meV, 
the MDC dispersion is almost vertical, with a weak S-like 
shape. We draw the attention to the fact that the hump 
shows a weakly positive dispersion close to the M point, 
with point of closest approach to the chemical potential 
at qy « O.SStt. This effect is due to the coupling of the 
(tt, 0) and (0, tt) points by self energy effects, and is a 
result of the fact that going towards F from the (tt, 0) 
point means going towards Y from the (0, tt) point at 
the (tt, 7r)-displaced wavevector. As a result of this, the 
weakening of the self energy effect along the M — F cut 
leads to a minimum in the hump dispersion at the M 
point. This effect was experimentally found. 

In Fig. ^1 we show our results for a cut parallel to 
the M — Y cut shown in Fig. keeping qy — 0.6tt con- 
stant. At low energies, the spectral evolution, seen on 
the left part of the figure, shows the typical BCS mixing 
between particle and hole states. Concentrating on the 
negative energy parts, again two branches are present, 
the peak branch and the hump branch, separated by a 
break in the EDC dispersion. Both branches now show 
considerable dispersion, but still overlap in momentum. 
The MDC dispersion changes from the low energy peak 
branch to the high energy hump branch at roughly the 
point where the intensity of the peak drops dramatically. 
Note that the EDC and MDC dispersions are consider- 
ably displaced relative to one another at high energies. 
Also at low energies, the MDC dispersion is stronger near 
the break region than the EDC dispersion. This effect in- 
creases when the residual width of the quasiparticle peak 
increases, and when convolution with the experimental 
resolution function is taken into account^ 

Finally, we discuss the cut along the nodal direction, 
shown in Fig. 1191 For this direction, the gap is zero as 
a consequence of d wave symmetry, and as a result the 
EDC dispersion should cross the Fermi energy. This is 
seen in the left panel of the figure. Note the very strong 
damping of the spectral peak as soon as it crosses the en- 
ergy region which corresponds to the break effect near the 
M point. Actually, the damping starts at slightly lower 
energies, due to the onset of node-node scattering pro- 
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FIG. 19: Left: Dispersion of the spectral intensity and line- 
shape in the nodal direction (F — y) as a function of momen- 
tum kx — ky = 0.257r...0.457r in steps of O.OItt from top to 
bottom. Right: the corresponding EDC (circles) and MDC 
(curve) dispersions. The kink is most clearly seen in the MDC 
dispersion. The low energy velocity is roughly half the high 
energy one. The high energy dispersion does not extrapolate 
to the Fermi surface crossing. 

cesses at an energy ^res, as can be seen in the left panel 
of Fig. ^1 The velocity renormalization for low energies 
and high energies differs by a factor of roughly two, both 
for EDCs and MDCs, in agreement with experiment. 
Finally, we also reproduce the experimental fact that the 
high energy dispersion does not extrapolate to the Fermi 
crossing. ^^-^^ Again, note some shift between the EDC 
and MDC dispersions at high energies due to the energy 
variation of the self energy. 

Clearly, the velocity break (kink) along the nodal di- 
rection and the break between the peak and hump (dip) 
near the M point are occurring in the same energy range 
between — fires — and — flres^ Em ■ This is an appeal- 
ing result of our theory, because it explains all features 
in the dispersion anomalies in the Brillouin zone seen by 
ARPES with a simple model. 

E. Tunneling spectra 

Knowing the spectral function, A(e,fc), throughout 
the zone, we are able to calculate the tunneling spec- 
tra given a tunneling matrix element Tg-. For simplicity, 
we present numerical results for the simple model, ne- 
glecting the continuum part of the spin fluctuation spec- 
trum. From the SIN tunneling current I{V), one obtains 
the differential conductance, dl/dV. As usual, we ne- 
glect the energy dependence of the SIN matrix element 
l^feP = '^eY.p\Ty\^AN{p,e), where An is the spectral 
function of the normal metal. The SIN tunneling current 
is then given by 

A{e,k){f{e)~f{e + eV)} 

k 

(38) 

We model the tunneling matrix element for two extreme 
cases: for incoherent tunneling we assume a constant 
|MgP = Mq, whereas for coherent tunneling we use 



FIG. 20: Differential tunneling conductance for SIN (left) 
and SIS (right) tunnel junctions for T = 40 K. Units are 
eMf for SIN and 2e'^T^ for SIS. Results for the coherent (full 
curves) and incoherent (dashed curves) tunneling limits are 
shown. The parameters are given in Table ITTl 

|MgP = lMf{cosk^ - coskyf^ Coherent tunneling 
in the c-axis direction is strongly enhanced for the M 
points in the Brillouin zone compared to the regions near 
the zone diagonal due to the matrix elements. ^'^ Our nu- 
merical results for SIN junctions are shown in Fig. I2(JI 
(left). In both cases, we observe a clear asymmetry, with 
a dip-hump structure on the negative bias side and a 
very weak feature on the positive side of the spectrum, 
as in experimentsi^2i2^ The low energy behavior of the 
tunneling spectrum in the coherent tunneling limit does 
not show the characteristic linear in energy behavior for 
d-wave, because the nodal electrons have suppressed tun- 
neling as a result of the matrix elements. The peak-dip- 
hump features, on the other hand, are not affected by the 
matrix elements, as they are dominated by the M point 
regions which are probed by both coherent and incoher- 
ent tunneling. 

For an SIS junction, the single particle tunneling cur- 
rent is given in terms of the spectral functions by 

I{V) = 2eY,\Ty\' r ^Aie,k)Aie + eV,p) 

J -OO ^'^ 

kp 

X {,/(£)-/(£ + el/)} (39) 

Again we show results for incoherent tunneling (iT'jTp'P = 
Tq) and for coherent tunneling with conserved parallel 
momentum, iTr-P — tsT? ( cos k^ — cos ky Y 5 r - .^"^ Our 
results are shown in the right panel of Fig. EOl 

All structures are symmetric around the chemical po- 
tential. The low energy part of the spectrum is strongly 
suppressed in the incoherent tunneling limit already, thus 
there is no big difference to the coherent tunneling limit 
there. At higher voltages, however, in the coherent tun- 
neling limit, we obtain negative differential conductance. 
Such an effect was observed recently in optimally doped 
Bi2Sr2CaCu208-5 break junctions."^^ We also observe 
negative behavior at higher bias in the coherent tunnel- 
ing limit, but note that in reality, the tunneling matrix 
element will have both coherent and incoherent contribu- 
tions (especially at higher voltages), and thus will be a 
weighted average of the dashed and full curves in Fig. 1201 
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In this case, most probably only the negative behavior 
below 100 meV will be observable. We note that the 
spin fluctuation continuum broadens the spectral func- 
tions and, as we show below, this leads to a positive re- 
sponse at higher voltages. 

We give approximate expressions for the SIS differen- 
tial conductances for zero temperature. In the incoherent 
limit, 



j{incoh) 



deNi£)N{e + eV) 



(40) 



-eV 



In the coherent tunneling limit, the tunneling matrix el- 
ement very effectively suppresses the nodal regions, thus 
only allowing for tunneling near the M point regions. In 
these regions, however, the dispersion is weak, so that 
we may approximate the spectral function by its value 
at the M point, Am{£)- Then, we obtain in the coherent 
tunneling limit 



j{coh)(y^ « _± / deAM{e)AM{e + eV) 

J~eV 



(41) 



with Ti — X^fcp'l^fcp'P- Note that two different quantities 
are probed in the two limits. In the incoherent limit, it is 
the density of states, and in the coherent limit, it is the 
spectral function at the M point of the Brillouin zone. 

It is easy to show by differentiating Eq.^Jthat the dif- 
ferential conductance can be negative, and furthermore, 
can approach a negative value for large voltages. The 
limiting behavior at high voltages in the incoherent tun- 
neling limit is proportional to A^(— )iV(-|-), where N(±.) 
is the density of states at large positive/negative ener- 
gies. If in the coherent tunneling limit the corresponding 
term proportional to Am{+)Am{—) is very small, then 
the main contribution comes from the region where either 
e « —Am or e + eV ~ Am, varying within a range of or- 
der Am around these values. It is easy to show that this 
contribution is negative. But as soon as incoherent con- 
tributions play any role, or if Am has a considerable in- 
coherent part, then their positive contributions will dom- 
inate at high voltages. Note that for SIN tunneling, the 
differential conductance is always positive definite. 



F. Self consistency issues 

When going towards underdoping, the spectral func- 
tion deviates considerably from the bare BCS spectrum. 
Self consistency issues become important then. 

Our studies have shown that the quasiparticle peak 
is always reasonably well separated in energy from the 
high energy incoherent part by a dip. By coupling elec- 
trons to the spin resonance mode, weight is shifted from 
the quasiparticle peak to the incoherent part which in- 
cludes the broad hump structure. Thus, when calculat- 
ing the self energy effects due to this coupling, only the 
quasiparticle peak part of the spectrum with its reduced 



weight will contribute to the sharp self energy features 
at energies affecting the quasiparticle peak. The inco- 
herent part of the fermionic spectrum, which is gapped 
by roughly the hump energy, will affect the low energy 
quasiparticle properties only in form of an effective high 
energy renormalization factor, which is constant up to 
energies comparable to the hump energy. This effective 
renormalization adds to the one due to coupling of elec- 
trons to the spin fluctuation continuum. Thus, we can 
concentrate on the renormalization equations following 
from the set of equations which includes the quasiparti- 
cle peak spectrum of reduced weight interacting with the 
spin fluctuation mode. In deriving these equations, we 
make use of the approximate equations for the renormal- 
ization functions derived above. 

The quasiparticle part of the Green's function has in 
this approximation at the M point the form 
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I^M + and Afj = (1 ± £,mIEm)/2. 

Here, Em is the measured peak position at the M point, 
and Fm is the quasiparticle peak width. The broadening 
of the off-diagonal spectra, 5m i is reduced compared to 
Fm due to d-wave symmetry. Using the approximative 
formulas from the last sections at e = —Em, we obtain 
(with a — q^Iq/i:) 
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Zm Z'^jflresEM ^res + 2Em 
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; Cm 
Ca/ = ^ - - 

^M 



Q.res + Em r , . 

Ui (46) 



Zm^tcsEm ^res + 2i?M 



where A^"^"* denotes renormalizations due to the spin fluc- 
tuation continuum and the incoherent part of the spec- 
tral function, and A^^^ the contribution coming from the 
nodal regions (these contributions are renormalizcd with 
the nodal renormalization factor Zn, which is smaller 
than Zm)- The last two equations merely express the 
measurable quantities Am and Cm as functions of the 
bare quantities Am and Cm- The first equation can be 
solved, giving for small fires and not too small g^Io a 



quasiparticle weight oc Qresi^res + 2-Bm)- Note that 

we derived this set of equations for the case where Fm 
is neglected, which describes the slightly underdoped re- 
gion. When fires becomes comparable to Fm , these equa- 
tions have to be modified. 

It should be remarked, though, that using these equa- 
tions in the absence of vertex corrections usually give 
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poorer results than those presented in this paper using 
bare Green's functionsi^ 



G. Bilayer splitting 

For bilayer compounds, the dispersion can be split into 
bonding (b) and antibonding (a) bands. Accordingly, the 
self energy for each band is defined as S^''^ and 
Similarly, the spin susceptibility is now a matrix in the 
bonding-antibonding indices, having elements diagonal 
iXaa, Xbb) and off-diagonal {xba, Xab) in the bonding- 
antibonding representation. The components of the spin 
susceptibility transforming even and odd with respect to 
the plane indices are given by Xe = Xaa + Xbb and Xo = 
Xab + Xba- For identical planes, we have Xaa = Xbb and 
Xab — Xba- The mcasurcd susceptibility is then given by 



X = xecos2^-Kxosm2^ (47) 

where d is the separation of the layers within a bilayer. 
If we write the self energy for a single layer as x * C 
(the hat denotes the 2x2 particle hole space), which is a 
functional of the spin susceptibility x and the Gor'kov- 
Green's function G, then we have formally for the two- 
layer system 



(48) 



For the resonance part, which only has a Xo component, 
this means that fermionic excitations of the antibonding 
band determine the self energy for the bonding band and 
vice versa. The calculations presented in this paper hold 
for the case of bilayer systems if we assume identical dis- 
persions for bonding and antibonding bands. Even small 
bilayer splittings of the order of 10 meV or less do not 
matter, as they do not qualitatively alter the spectral 
form of the self energy. For larger bilayer splittings, the 
self energy is larger for the bonding band, because it is 
determined by the van Hove singularity near the chemical 
potential in the antibonding band. Thus, stronger renor- 
malizations are expected in the bonding band for this 
case, which tends to decrease the bonding-antibonding 
splitting. This effect of reducing the bilayer splitting 
should be strongest in underdoped compounds, where the 
effect of the resonance mode is strongest. In overdoped 
compounds, the bilayer splitting should be less affected 
by spin fluctuations. Our prediction is that if a bilayer 
splitting is observed, then the peak-dip-hump structure 
should be stronger for the bonding band with the higher 
binding energy peak. This is consistent with the data of 
Ref. ^5. The onset of strong fermionic damping should 
be independent of the band index, as it is given by scat- 
tering to the nodes, and thus occurs at the fixed energy 

^res ■ 

In this paper, we have elected not to explicitly include 
bilayer splitting effects in our calculations. The primary 



reason is that although all ARPES groups now detect 
the presence of bilayer splitting for heavily overdoped 
samples, the various groups disagree on its presence for 
optimal and underdoped samples—. Recently, we have 
performed calculations including bilayer splitting and are 
able to reproduce a number of unusual spectral anoma- 
lies seen in heavily overdoped ARPES spectraSi. These 
calculations further confirm the picture advocated in this 
paper, in that the spectral anomalies imply a mode which 
has odd symmetry with respect to the layer index of the 
bilayer, a unique property of the magnetic resonance ob- 
served by neutrons. For further details, the reader is 
referred to Ref. 1671 



H. Doping dependence 

In this section, we deal with the doping dependence of 
the spectral lineshape near the M point of the Brillouin 
zone. As there are many parameters which change with 
doping in different ways, it could turn out to be a mean- 
ingless task to adjust all of those parameters and at the 
same time make a sensible prediction. But, fortunately, 
all changes with doping lead to spectral changes which go 
in the same direction. This 'fortuitous' accident allows 
us to make some general predictions from the theory we 
use. To see this, we turn again to Figs. IlllandlT^ From 
there we see that the quasiparticle weight decreases with 
decreasing fires /'^m, and with increasing coupling con- 
stant g^WQ. The quasiparticle scattering rate increases 
with decreasing fi^es/AM- And the hump energy dis- 
perses to higher binding energies for increasing coupling 
constant and increasing ^m- Thus, in our model, going 
from overdoping to underdoping amounts to a decreasing 
quasiparticle weight, an increasing quasiparticle scatter- 
ing rate, and an increasing hump binding energy. 

The important parameter, as we see from this study, 
is the ratio fires /^m, the ratio of the mode energy to 
the maximal superconducting d-wave gap. We distinct 
two regions: the first, where fires/ > 1, and the sec- 
ond, where fires/ < 1- The situation is schematically 
shown in the phase diagram in Fig. 1211 The curves shown 
are calculated using the formulas (we relate Tc to the hole 
doping level in the CU-O2 planes in the usual manner^) 



Tc = 95 K(l -82. 6(p- 0.16)2) 
Am = 38meV(l-9.1(p-0.16)) 
fires = 40 meV(l-82.6(p- 0.16)2) 



(49) 
(50) 
(51) 



All these quantities approach zero on the overdoped side 
sX p — 0.27. Optimal doping corresponds to p = 0.16. 
Note that fires = 4.9Tc in agreement with Ref. IstI The 
Am variation was based on ARPES datai^Mi 

As can be seen, the separation between overdoped and 
underdoped regions roughly coincides with the regions 
where fires > and fires < Am, respectively. The dip 

A4. As A 4 is about the same 



onset is given by fires + ' 
as Am, we have shown in Fig. [2] the line for fir 
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FIG. 21: In the dark gray shaded region, corresponding to 
overdoping, quasiparticles peaks are well defined. In the light 
gray shaded region, corresponding to underdoping, the peak 
weight is strongly reduced, and an incoherent part due to scat- 
tering from the spin fluctuation mode is dominant. The reso- 
nance energy, shown as a thick line, is bounded from above by 
twice the maximal gap energy, Qres < '2 Am, and approaches 
it on the overdoped side. The position of optimal doping, at 
maximal Tc and Qres , roughly coincides with the point where 
Am as a function of doping crosses fires- 



as a dashed line, which determines the position of the 
dip fairly accurately. The continuum in the spin fluctu- 
ation spectrum only affects electrons above 2A/,,, which 
is near or above the dip energy. One important observa- 
tion is that the point of optimal doping for a CU-O2 plane 
roughly corresponds to the point where fires /^m — 1- 
Thus, region I of Fig. is relevant to overdoped mate- 
rials, and region II to underdoped materials. Another 
experimental observation is that this ratio never exceeds 
the value two. This is expected for an excitonic collective 
mode below a continuum edge^ 

For a quantitative theory of the doping dependence the 
self consistency issue becomes important. The coherent 
quasiparticle weight and the quasiparticle linewidth are 
given by. 
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where Zm = Zm{—Em) is the only quantity not avail- 
able from experiment. We can eliminate it, to obtain the 
relation 
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Note that experimentally both zmEm and (possibly) 
-Da scale with ksTcMiSl. So, the quasiparticle width is 



dominated by the difference Em — ^res- Quasiparticles 
are sharp at the overdoped side where Em < f^resj and 
an onset of quasiparticle scattering as a function of un- 
derdoping takes place when Em — ^res- This point is 
slightly beyond optimal doping. 



IV. COUPLING TO THE ZONE BOUNDARY 
HALF BREATHING OPTICAL PHONON 

The sharp structure in the dispersion needs an expla- 
nation in terms of an almost dispersionless feature which 
couples to the electrons. Numerous phonons modes are 
seen in inelastic neutron scattering in high Tc cuprates. 
Most of them do not show indications of strong coupling 
to electrons. Two special types of phonons have attracted 
attention: the Cu-0 buckling mode, which is attractive in 
the d-wave channel)^i^^i''°i'^^i''^i'''^ and the Cu-0 breathing 
mode, which is repulsive in the d-wave channel iSi£2iiiiS 
Typically, the absolute values of the pairing interactions 
in the Big ('d-wave') channel for both types of vibra- 
tions are smaller than 0.1 eV, in the Aig ('s-wave') 
channel about 0.5-1 eV; for spin fluctuations, the cor- 
responding numbers are in the d-wave channel 0.5-1 eV 
and in the s-wave channel 1-2 eV.?3 The total electron- 
phonon coupling constant in the s-wave channel amounts 
to \s w 0.4 - 0.6/3,74,75,76,77 g^^^^ ^j^g d-wave channel 

to Ad ~ 0.3i2Si£i Thus, phonons are not likely to be re- 
sponsible for the high transition temperature. 

It was argued recently that strong coupling of elec- 
trons to the zone boundary half breathing phonon may 
be responsible for the anomalies in the dispersion. It is 
known for some time that this phonon shows a disper- 
sion which is strongly renormalized midway between the 
zone boundary and the zone center when entering the su- 
perconducting state. These findings show that the zone 
boundary half breathing phonon is affected by super- 
conductivity. It was suggested to be responsible for the 
renormalizations of the dispersion observed in ARPES.^— 
This zone boundary half breathing longitudinal optical 
phonon is a Cu-0 bond stretching mode with an energy 
between 50 and 100 meV. Its dispersion is very strong 
in the middle of the branch, and it was suggested that a 
discontinuity develops there in the metallic statei^^ The 
first measurements concentrated on lanthanum cuprates, 
but recently YBa2Cu307_a was also studiedi ^^i^° The 
displacements involve oscillations of the oxygen atoms 
in phase between the two planes in the bilayer. The 
results for optimally doped YBa2Cu307_5 are the fol- 
lowing: The dispersion of the zone edge mode in the 
superconducting state shows a 'break' at (0, 7r/2) (and 
equivalent points), with an almost dispersionless branch 
(at ~ 55 meV) between q = (0, 7r/2) and q = (0,7r), 
and a dispersive branch (68 meV to 72 meV) between 
g = (0,0) and q = (0, 7r/2)i7?i7?i?? Experimental inves- 
tigation showed that the dispersionless branch extends 
over a region tt/2 < q^ < tt, — O.Itt < q^ < O.Itt (and 
analogously for q^ and qy interchanged) i2S The disper- 
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tered by these phonons due to the Pauh exclusion princi- 
ple. The important points here are near the nodal regions 
for the electrons. Thus, the strongest effects are expected 
near (but not necessarily at) the nodes, not near the Ap- 
points, in contrast to what experiment shows. 

It is possible that both processes play a role and dom- 
inate in different regions of the Brillouin zone. Phonons 
would then play some role for nodal electrons. 



CONCLUSIONS 



FIG. 22: As in Figs. Eld) and0] but the white areas now 
denote the experimentally determined regions where the dis- 
persionless half breathing zone boundary phonon mode is 
present. For nodal wave vectors (left) the main self energy 
contribution comes from node-node scattering processes at 
small energies (near —flphon)- At the M point (right) the 
self energy effects are negligible due to geometric restrictions. 
Only higher anharmonic terms (with two phonon processes 
which add up to a (tt, tt) wavevector) could contribute at 
fermionic wavevectors near the M point. 



sionless branch was only observed for bond stretching 
along the a direction (perpendicular to the chains). The 
dispersions of the longitudinal bond-stretching phonon 
branches were found to show no apparent temperature 
dependence^ The phonon intensity was found to show 
significant temperature dependence below Tc7^ Phonon 
weight is transferred from a position halfway to the zone 
boundary (in a range between 55 meV and 70 meV) to 
the zone center (70-75 meV) and the zone boundary (50- 
55 meV). This transfer sets in at Tc and increases with 
decreasing temperature. 

The coupling strength, gbio), goes to zero for small 
momentum transfer q. Furthermore, in the model of 
Ref. Isil the coupling vanishes near the q = (tt, tt) point, 
thus having minimum impact on the electrons near the 
M point of their Brillouin zone. This is in stark contrast 
to the resonance mode model, and can certainly not ex- 
plain the effects at the M points. It is, however, possible 
that they contribute to the renormalization of the nodal 
dispersion. The maximal coupling strength was theoret- 
ically estimated to gb ~ 0.04eVp2,, but in some models is 
enhanced by vertex corrections^ 

In Fig. 1221 we point out an important difference to the 
magnetic mode. The magnetic mode is peaked around 
(tt, tt), whereas the spectral density of half breathing 
phonon is peaked around the points (0,7r) and (7r,0). 
Because for the imaginary part of the self energy only ex- 
citations near the Fermi surface are important, there are 
geometric restrictions for the possible scattering events. 
In the case of the magnetic mode, scattering was domi- 
nated by processes connecting the M points of the Bril- 
louin zone, and these scattering processes are enhanced 
by the presence of a van Hove singularity close to the 
chemical potential. In the case of the half breathing 
phonon, the M point electrons are very ineffectively scat- 



We suggest that the van Hove singularity at the M 
point of the Brillouin zone plays an important role in de- 
termining the self energy effects observed in ARPES and 
tunneling experiments. The picture can be understood 
as follows: the quasiparticle dispersion is fairly flat near 
the M points of the Brillouin zone, with a large effective 
mass in the M — P direction, and is close to the chemical 
potential. Because the continuum part of the spin fluc- 
tuation spectrum is gapped up to energies ~ 70 meV, 
the scattering at low energies is dominated by scattering 
processes accompanied by emission of a spin fluctuation 
mode (which lies below the gapped continuum at ~ 40 
meV and has a sharp energy). This coupling leads to 
cusps in the energy dependence of the self energy in the 
range of 70-80 meV due to the effect of the van Hove sin- 
gularities at the M and A points. Because of the finite 
width of the spin fluctuation mode in momentum, there 
are traces of these cusps for electrons at all positions near 
the Fermi surface. The position in energy of these cusps 
are determined by electrons near the M and A points 
only, which explains the isotropy around the Fermi sur- 
face of the energy scale of 70-80 meV where kink fea- 
tures in the dispersion are observed. The intensity of 
this self energy effect is determined by the intensities of 
the spin fluctuation mode at such momenta q which con- 
nect the electron with momentum k to electrons near the 
M point region. Thus, q — {k — modulo (G), where 
G is a reciprocal lattice vector. This intensity is large 
for k « fcM, but smaller for k ~ k^- This explains the 
strong anisotropy of the magnitude of the effect around 
the Fermi surface. 
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Note added in proof. After submitting this paper, we 
became aware of an experimental paper which claims 
no momentum anisotropy in the linewidth for overdoped 



22 



compounds^. This result is actually consistent with our 
picture, in that for overdoped compounds, the spectral 
peak lies inside the scattering rate gap, which can be 
appreciated from Fig. 21. 



where D^{e — ien) and G^{€ — ium) are smooth func- 
tions (except at ujm = 0, which is treated separately, see 
below). So, the self energy Ea. lAll has the two alternative 
equivalent forms (the first form was found in Ref. [s^, 



APPENDIX A: SELF ENERGIES 

In this appendix, we derive self energy expressions 
which allow us to evaluate the real part of the self en- 
ergy analytically in several special cases, and have the 
numerical advantage of having eliminated all principal 
value integrals. The procedure is a generalization of a 
method developed by F. Marsiglio et al^ The self en- 
ergy is given by. 



(Al) 



where D — — In equilibrium, the Keldysh components 
are given by the simple expressions 
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where B^^^q and A^^k are the bosonic and fermionic spec- 
tral functions, and b^, their corresponding distribution 
functions, respectively. Note, that the Keldysh compo- 
nents are purely imaginary. 

In evaluating these integrals, the only numerical com- 
plication comes from the real parts of S due to principal 
value integrals. We can eliminate those by using the fol- 
lowing trick: Note that in equilibrium the identities 
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hold which are easy to check. The convenient feature is 
that the second lines in Eas. lA4l and lA5l can be converted 
into Matsubara sums by noting that D^Gf_^ is an ana- 
lytic function in the lower lj half plane, and analogously 
D^Gf_^^ analytic in the upper half plane. Thus, 
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where the population factor p^ ,:_^ is given by, 
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Note that the terms containing Matsubara sums are pure 
real quantities. 

Let us examine the simple case 



= 2n6ie - a) 
which gives, using the second expression 
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Finally, for the case that the bosonic mode has the simple 
form 
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the first expression leads to 
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The last sum can be performed for the case of a Green's 
function of the form 
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